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SUMMARY 
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A  basic  language  computer  program  is  presented  which  calculates  the  mois¬ 
ture  absorption  of  epoxy  resin  and  fibre  composite  sheets  from  the  sheet  thick¬ 
ness,  the  diffusivity  and  the  relative  humidity  of  the  environment.  One¬ 
dimensional  Fickian  diffusion  is  assumed  in  the  calculations.  The  program  gives 
the  water  absorbed  by  the  material  in  a  given  exposure  time,  together  with  the 
times  for  the  material  to  absorb  95%  and  99%  of  the  maximum  water  uptake.  In 
addition,  the  distribution  of  water  through  the  laminate  thickness  is  calculated 
and  displayed  graphically.  The  programme  can  be  used  for  any  resin  or  fibre 
composite  in  which  the  primary  moisture  diffusion  mechanism  is  bulk  diffusion 
through  the  resin.  Examples  for  typical  fibre  reinforced  epoxy  composite  lami¬ 
nates  are  given  and  the  implications  are  discussed.  Diffusion  data  for  three 

* 

materials  at  room  temperature  and  two  materials  at  50  C  are  listed  for  use  in 
the  program. 
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1  INTRODUCTION 

Fibre  reinforced  plastics,  because  of  their  high  strength  and  stiffness  combined 
with  their  low  density  and  the  low  cost  of  the  finished  components,  are  being  used 
increasingly  in  airframe  structures.  In  a  service  environment,  fibre  composite  struc¬ 
tures  will  be  exposed  to  ultra-violet  radiation,  lightning  strikes,  acoustic  fatigue, 
corrosive  fluids  and  atmospheric  moisture  and  consideration  must  be  given  to  their  effect 
on  mechanical  properties.  In  this  Memorandum  only  absorbed  moisture  is  considered. 

The  current  generation  of  epoxy  resins  used  in  high  performance  fibre  composites 

absorb  water  from  the  atmosphere.  The  immediate  effect  of  this  is  a  swelling  of  the 

resin  which  counteracts  to  some  extent  the  shrinkage  during  the  curing  process  and  it  can 

result  in  significantly  reduced  residual  thermal  strains  in  laminates1.  Water  absorption 

by  the  epoxy  resins  also  leads  to  a  reduction  in  the  glass  transition  temperature  and  to 

a  softening  of  the  resin  with  a  loss  of  resin  stiffness  and  strength,  particularly  at 
2 

elevated  temperatures  .  In  the  longer  term  the  resin  may  be  permanently  damaged  and  a 

3 

further  degradation  m  properties  may  result  .  These  degraded  resin  properties  manifest 
themselves  in  the  fibre  composite  as  a  loss  of  performance  in  the  resin  dominated  proper¬ 
ties  such  as  reductions  in  strength  and  stiffness  under  shear  loading,  compressive  load- 

2  4  5 

ing  and  loading  perpendicular  to  the  fibres  *  *  .  In  some  composites,  such  as  those 
reinforced  with  polymer  fibres  or  glass  fibres,  the  fibre  properties  may  also  be  degraded 
by  moisture  absorption. 

The  quantity  of  water  absorbed  by  a  laminate  is  thus  of  considerable  importance,  in 
particular  to  the  designer  when  setting  design  limits  for  structures  operating  in  moist 
environments.  In  this  Memorandum  a  basic  computer  program  is  presented  for  the  calcula¬ 
tion  of  the  quantity  of  water  absorbed  by  a  fibre  composite  or  resin  sheet  for  a  given 
exposure  time  in  an  environment  of  known  relative  humidity.  It  is  assumed  that  the  rate 
controlling  mechanism  is  one-dimensional  Fickian  diffusion  through  the  resin,  and  that 
diffusion  of  moisture  in  the  fibres  or  along  the  fibre-resin  interface  is  negligible. 

The  equilibrium  quantity  of  water  absorbed  in  the  environment  is  calculated  as  well  as 
that  in  the  specified  time.  In  addition  the  time  taken  for  the  laminate  to  absorb  95% 
and  99%  of  the  maximum  water  uptake  in  the  environment  is  calculated.  The  program  also 
permits  a  profile  of  the  distribution  of  water  through  the  laminate  thickness  to  be 
plotted  for  the  given  conditions.  The  only  input  requirements  are  the  laminate  thick¬ 
ness,  exposure  time,  relative  humidity  of  the  environment,  the  saturation  water  content 
(in  a  100%  RH  environment)  the  diffusivity  of  water  in  the  material  and  the  initial 
moisture  content.  Dif f usivities  of  three  materials  at  two  temperatures  are  given  in 
section  5,  together  with  the  saturation  water  contents  of  the  materials.  The  program 
was  written  specifically  to  run  on  a  CBM  PET  computer  with  a  usable  RAM  (random  access 
memory)  of  at  least  8K  bytes.  However,  the  program  could  easily  be  modified  to  run  with 
any  basic  language  compiler. 

2  MOISTURE  ABSORPTION  THEORY 

In  a  glass,  carbon  or  boron  fibre  reinforced  epoxy  resin,  moisture  is  absorbed  by 

2  6 

the  resin;  the  fibres  do  not  absorb  moisture.  Most  of  the  evidence  in  the  literature  * 
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suggests  that  moisture  is  absorbed  by  a  bulk  diffusion  mechanism  in  the  resin  and  that 
for  sheet  laminates  the  rate  of  absorption  3c/3t  can  be  described  by  Fick's  second 
law  of  diffusion  in  one  dimension  ( ie  through  the  thickness): 


Sc 

3t 


(1) 


where  c  is  the  moisture  concentration  at  time  t  and  x  is  the  distance  from  the 

laminate  surface.  D  is  the  diffusion  coefficient  (dif fusivity)  in  the  x  direction; 
x 

this  is  a  thermally  activated  parameter  of  the  Arrhenius  type  and  varies  with  temperature 

such  that  diffusion  occurs  more  rapidly  at  higher  temperatures.  The  fit  between  Fickian 

diffusion  theory  and  experiment  is  generally  good  but  not  perfect  and  there  is  evidence 
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that  in  some  circumstances  there  can  be  large  discrepancies  9  .  These  may  be  due  to 
anomalous  diffusion  processes  such  as  more  rapid  diffusion  along  the  fibre-resin  inter¬ 
faces  (or  fibres  in  the  case  of  polymer  fibre  reinforced  resins)  which  is  not  included 
in  bulk  diffusion  theory.  In  addition,  the  diffusivity  may  also  be  a  function  of  lami¬ 
nate  thickness,  fibre  orientation,  internal  stress  or  moisture  concentration,  in  which 
case  equation  (1)  would  be  inapplicable  and  the  problem  would  be  significantly  more 
complex.  At  present,  however,  Fickian  diffusion  theory  provides  a  simple  technique  for 
describing  the  approximate  moisture  absorption  characteristics  of  resins  and  fibre  com¬ 
posite  laminates  and  this  is  the  basis  of  the  calculations  in  this  work. 

The  solution  of  equation  (1)  has  been  given  by  Shen  and  Springer^  as: 


on 

j=0 


2j  +  1 


sin(2j  +  1)ttx 


exp 


(2 j  +  1)2w2D  t 
x 


(2) 


where  is  the  initial  concentration  at  zero  time  and  c  is  the  maximum  concentra- 

U  m 

tion  in  the  given  environment  (both  per  unit  volume).  Equation  (2)  thus  gives  the  mois¬ 
ture  concentration  through  the  laminate  thickness  h  at  any  time  t  and  position  x  . 
This  equation  has  been  used  in  this  work  to  calculate  moisture  profiles  across  laminate 
thicknesses.  The  summation  was  performed  for  up  to  j  -  100  ,  but  the  series  usually 
converges  for  j  <  10  . 

The  total  moisture  content  per  unit  volume  in  a  laminate  is  determined  by  integrat 
ing  equation  (2)  over  the  laminate  thickness 


0 


(3) 


Shen  and  Springer  have  given  an  approximate  solution  to  equation  (3)  which  fits  the 
appropriate  boundary  conditions^: 
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(4) 


where  m  is  the  total  moisture  content  per  unit  volume  in  the  laminate  at  time  t  , 

is  the  maximum  moisture  content  per  unit  volume  in  the  given  environment  and  m^  the 
moisture  content  per  unit  volume  at  zero  time.  In  practical  applications  it  is  the  per¬ 
centage  moisture  content  M  that  is  of  interest,  as  given  by 


M 


(5) 


where  W  is  the  mass  of  moist  material  and  that  of  the  dry  material.  Since 

W  *  +  m  ,  it  is  apparent  that: 
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m  -  m 
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m 
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M  -  M. 
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(6) 


where  is  the  initial  and  the  maximum  percentage  moisture  content.  Current 

evidence  on  typical  epoxy  based  laminates  supports  the  idea  that,  although  is  a 

strong  function  of  temperature,  is  independent  of  temperature  (at  least  within  the 

normal  range  of  working  temperatures) .  Equations  (4)  and  (6)  have  been  used  in  this  work 
to  calculate  the  total  moisture  content  M  in  the  specified  environment  for  a  known  time 
t  and  laminate  thickness  h  .  For  large  t  ,  the  quantity  G  in  equation  (4)  tends  to 
unity,  only  reaching  this  at  infinite  time.  To  make  an  estimate  of  the  time  for  M  to 
approach  it  is  thus  necessary  to  set  G  to  less  than  unity.  In  this  work  the  time 

for  M  to  reach  95%  and  99%  of  was  calculated  by  the  appropriate  substitution  in 

equation  (4). 

The  maximum  moisture  content  in  an  environment  with  <{>%  relative  humidity  is  given 
by  the  empirical  relationship^: 

M  =  a($/100)b  (7) 

m 


where  a  and  b  are  experimentally  determined  constants  for  each  material  (a  is  the 

saturation  moisture  content  in  %  of  the  material  in  a  100%  RH  environment).  For  many 

materials,  b  in  equation  (7)  can  be  set  approximately  to  unity  and  only  the  value  of 

9  10 

the  parameter  a  varies  with  material  type  ’  .  This  approximation  usually  provides  an 

adequate  description  of  moisture  content  as  a  function  of  relative  humidity,  although 

deviations  frequently  occur  in  very  wet  environments^.  Thus  in  this  work  b  was  set 

9  1  1 

to  unity  and  values  of  the  parameter  a  were  obtained  from  various  sources  * 


3  INPUT  DATA  REQUIREMENTS 

Before  the  program  may  be  executed,  the  following  parameters  must  be  known: 

(a)  The  laminate  thickness  in  mm. 

(b)  The  relative  humidity  of  the  environment  in  %  . 
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(c)  The  exposure  time  in  days. 

(d)  The  initial  moisture  content  in  %  . 

(e)  The  name  of  the  resin  for  reference  purposes. 

.  .  2  -1 

(f)  The  diffusivity  of  the  material  at  the  required  temperature  in  mm  s 

(g)  The  saturation  water  content  of  the  material  in  %  (in  a  100%  RH  environment). 

4  THE  COMPUTER  PROGRAM 

A  complete  listing  of  the  program  is  given  in  Table  I.  This  should  be  entered  into 
the  computer1 s  memory ,  then  the  command  RUN  followed  by  RETURN  will  cause  the  program  to 
be  executed. 

Between  lines  180  and  260  of  the  program  the  operator  is  requested  to  enter  the 
name,  diffusivity  and  maximum  moisture  uptake  of  the  resin  based  material.  At  line  190 
the  relative  humidity  of  the  environment  should  be  entered  on  request.  In  addition  the 
exposure  time,  laminate  thickness  and  initial  moisture  content  should  be  entered  as 
directed.  Equation  (4)  is  used  to  calculate  the  moisture  content  M  in  the  specified 
conditions;  this  appears  as  lines  280-290  in  the  program.  Lines  300-330  are  concerned 
with  printing  the  results  of  the  initial  calculations  on  the  vdu  screen.  Lines  340-350 
and  410-420  make  use  of  equation  (4)  to  calculate  the  time  taken  for  the  material  to 
reach  95%  and  99%  respectively  of  the  maximum  moisture  content  .  This  information 
is  printed  out  through  lines  370,  400,  450  and  470.  At  line  490  the  operator  is  asked 
whether  a  moisture  profile  across  the  laminate  thickness  is  required.  The  GET  command 
in  line  500  responds  to  a  positive  reply  by  initiating  the  graph  plotting  subroutine 
from  line  570  onwards.  Any  other  response  leads  to  termination  of  the  programme  at 
line  550,  but  a  re-run  with  the  same  material  and  environmental  conditions  can  be 
triggered  by  a  positive  response  to  the  GET  statement  of  line  530.  Output  on  either  the 
computer  vdu  or  an  external  printer  can  be  selected  by  the  appropriate  command  at  lines 
580-590. 

In  lines  670-770  the  distribution  of  moisture  across  the  laminate  thickness  is 
calculated  using  equation  (2).  Thirty  values  are  calculated  at  equal  spacings  across 
the  thickness.  The  rest  of  the  programme  is  a  graph  plotting  routine  which  displays  the 
water  content  of  the  laminate  versus  the  distance  across  the  laminate  thickness. 

Automatic  scaling  and  labelling  of  the  axes  is  incorporated  in  the  routine.  When  the 
plot  is  complete  a  RETURN  command  at  line  2110  returns  control  to  line  510  and  at  line 
520  the  operator  can  request  a  re-run  with  the  same  material  and  environmental  conditions 
by  typing  Y  .  Any  other  response  leads  to  the  termination  of  the  program. 

5  TYPICAL  RESULTS 

Dif fusivities  and  values  of  the  parameter  a  for  three  epoxy  resin  based  mater¬ 
ials  at  room  temperature  and  two  at  50°C  are  listed  in  Table  2  (fibre  volume  fractions 
are  approximately  0.6  for  all  five  cases).  This  information  may  be  loaded  into  the 
program  as  directed  (see  section  4  for  details). 

The  following  data  for  an  XA-S/914  laminate  exposed  at  room  temperature  was 
entered  as  described  in  section  4.  Use  was  made  of  the  data  in  Table  2. 
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5 . 1  Example  1 

(a)  Laminate  thickness  *  2  nan. 

(b)  The  relative  humidity  of  the  environment  =  65%. 

(c)  Exposure  time  =  30  days. 

(d)  Initial  moisture  content  =  0%. 

(e)  Resin  name  =  914. 

The  output  for  this  material  is  summarised  below: 

MAXIMUM  MOISTURE  CONTENT  -  1.125% 

MOISTURE  CONTENT  =  0.6653%  AT  30  DAYS 

TIME  TO  REACH  95%  OF  MAXIMUM  WEIGHT  INCREASE  IS  150.19  DAYS 
TIME  TO  REACH  99%  OF  MAXIMUM  WEIGHT  INCREASE  IS  266.47  DAYS 

The  moisture  profile  across  the  laminate  thickness  is  given  in  Fig  1. 

5.2  Example  2 

(a)  Laminate  thickness  -  10  mm. 

(b)  The  relative  humidity  of  the  environment  -  65%. 

(c)  Exposure  time  =  100  days. 

(d)  Initial  moisture  content  =  0%. 

(e)  Resin  name  =  914. 

The  output  for  this  material  is  summarised  below: 

MAXIMUM  MOISTURE  CONTENT  =  1.125% 

MOISTURE  CONTENT  =  0.2016%  AT  100  DAYS 

TIME  TO  REACH  95%  OF  MAXIMUM  WEIGHT  INCREASE  IS  10  YEARS  AND  104.9  DAYS 
TIME  TO  REACH  99%  OF  MAXIMUM  WEIGHT  INCREASE  IS  18  YEARS  AND  91.75  DAYS 

The  moisture  profile  across  the  laminate  thickness  is  given  in  Fig  2. 

6  DISCUSSION 

The  computer  program  presented  in  section  4  permits  the  moisture  absorption  behav¬ 
iour  of  thin  sheets  of  resins  or  fibre  reinforced  resins  to  be  calculated.  Example  data 
is  included  in  section  5  which  permits  the  behaviour  of  three  materials  at  room  tempera¬ 
ture  and  two  materials  at  50°C  to  be  calculated,  although  the  moisture  absorption 
behaviour  may  be  calculated  for  any  sheet  material  if  the  diffusivity  at  the  required 
temperature  and  the  saturation  moisture  content  are  known. 

In  a  service  environment,  both  the  relative  humidity  and  temperature  of  the 

environment  will  be  changing  daily  and  seasonally,  so  that  both  the  diffusivity  and  the 

maximum  moisture  content  will  vary.  The  relative  humidity  may  vary  from  dry  to  wet, 

0  to  100%  RH,  although  in  a  typical  North  European  environment  this  range  would  probably 

be  smaller.  The  operating  temperature  of  a  typical  military  aircraft  may  vary  from 

-50  to  +I20°C.  Thus  the  maximum  moisture  content  of  a  resin-based  structural  component 

—8 

will  vary  considerably  and  the  diffusivity  of  the  material  also,  from  less  than  10  to 
-7  2-1 

greater  than  10  /mm  s  .  However,  an  average  quasi-steady  state  moisture  content  is 
reached  after  long  exposure  times,  typically  about  10  years  for  panels  more  than  a  few  ntn 
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thick,  with  only  minor  fluctuations  superimposed  .  The  moisture  distribution  in  the 

surface  layers,  however,  is  more  changeable  because  they  can  respond  more  rapidly. 

Although  the  computer  program  is  only  applicable  for  constant  relative  humidity  environ¬ 
ments,  this  behaviour  is  demonstrated  in  the  example  in  section  5.2  for  a  typical  North 
European  environment  of  65%  relative  humidity.  In  this  example  a  10mm  thick  CFRP  lamin¬ 
ate  was  shown  to  take  over  10  years  to  absorb  95%  of  the  maximum  moisture  content,  but 
after  100  days  continuous  exposure  only  the  surface  layers  had  absorbed  significant 
quantities  of  moisture,  the  interior  of  the  laminate  still  being  dry. 

7  CONCLUSIONS 

A  basic  computer  program  was  presented  and  discussed  that  allows  the  moisture  con¬ 
tent  to  be  calculated  in  a  thin  sheet  of  resin  based  material  exposed  to  an  environment 
of  known  temperature  and  relative  humidity,  for  a  specified  time.  The  program  also  per¬ 
mits  the  calculation  of  the  time  taken  for  the  material  to  absorb  95%  and  99%  of  the 
maximum  moisture  content  in  the  specified  environment.  The  distribution  of  moisture 
through  the  sheet  thickness  may  also  be  plotted.  The  programme  can  be  used  for  any 
resin  or  fibre  composite  in  which  one-dimensional  Fickian  diffusion  adequately  describes 
the  moisture  absorption  mechanism  through  the  sheet  thickness.  Examples  for  typical 
fibre  composites  are  given.  Fibre  composite  laminates  more  than  a  few  mm  thick  were  shown 
to  take  many  years  to  reach  saturation  and  in  the  short  term  only  the  surface  layers 
responded  rapidly  to  the  exposure  environment. 
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Table  1  (continued) 
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Tab 1 e  I  ( c  on  t i n ued ) 


Table  1  (continued) 


Table  2 


RESIN  DIFFUSrVITY  AND  MOISTURE  CONTENT  DATA 


Temperature 

Fibre  type  and 

Resin  type  and 

Dif fusivity 

min^  S“1 

degrees  C 

manufacturer 

manufacturer 

_ _ 

23 

XA-S 

Courtaulds 

BSL-914C 

Ciba-Geigy 

9.4  x  I 0~8 

1 

■  23 

XA-S 

Courtaulds 

Code  69 

Fothergill  and 

5.04  x  io~8 

Harvey 

t 

i 

I  23 

A-S 

Courtaulds 

3501-5 

Hercules 

o 

O' 

X 

o 

1 

50 

l 

HT-S 

3501-5 

1  .2  X  10-7 

Courtaulds 

Hercules 

50  ! 

XA-S 

Courtaulds 

BSL  914C 
Ciba-Geigy 

1.26  x  io'7 

percent  > 


Source 


Ref 

9 

Ref 

9 

Ref 

1 1 

Ref 

12 

Ref 

12 

(*  estimated  from  known  moisture  content  in  a  65%  RH  environment) 
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10 
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Fig  1 


Fig  1  A  typical  output  for  a  2m  thick  XAS/914  laminate  exposed 
for  30  days  at  65%  RH  23°C 
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